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ABSTRACT 



G. P. Smith,' K. Umetsu,^ T. Futamase/ 



We present a joint weak-lensing/X-ray study of galaxy cluster mass-observable scaling relations, 
motivated by the critical importance of accurate calibration of mass proxies for future X-ray missions, 
including eROSITA. We use a sample of 12 clusters at z ~ 0.2 that we have observed with Subaru and 
XMM-Newton to construct relationships between the weak-lensing mass (M) and three X-ray observ- 
ables: gas temperature (T), gas mass (Mgas), and quasi-integrated gas pressure (Ix) at overdensities 
of A = 2500, 1000, and 500 with respect to the critical density. We find that Mgas at A < 1000 
appears to be the most promising mass proxy of the three because it has the lowest intrinsic scatter 
in mass at fixed observable, ainM — 0.1, independent of the cluster dynamical state. The scatter in 
mass at fixed T and Ix is a factor of '^ 2 — 3 larger than at fixed Mgas, which are indicative of the 
structural segregation that we find in the M — T and M — Ix relationships. Undisturbed clusters are 
found to be ~ 40% and ^ 20% more massive than disturbed clusters at fixed T and Ix respectively 
at '^ 2(T significance. In particular, A 1914 - a well-known merging cluster - significantly increases the 
scatter and lowers the normalization of the relation for disturbed clusters. We also investigated the 
covariance between intrinsic scatter in M — Afgas and M — T relations, finding that they are positively 
correlated. This contradicts the adaptive mesh refinement simulations that motivated the idea that 
Ix rnay be a low scatter mass proxy, and agrees with more recent smoothed particle hydrodynamic 
simulations based on the Millennium Simulation. We also propose a method to identify a robust mass 
proxy based on principal component analysis. The statistical precision of our results is limited by the 
small sample size and the presence of the extreme merging cluster in our sample. We therefore look 
forward to studying a larger, more complete sample in the future. 

Subject headings: Cosmology: observations - dark matter - galaxies: clusters: general - gravitational 
lensing: weak - X-rays: galaxies: clusters. 



1. INTRODUCTION 

Galaxy clusters are the largest virialized objects in 
the universe; they formed from high amplitude peaks 
of the primordial density field. Clusters therefore oc- 
cupy the high mass exponential tail of the dark mat- 
ter halo mass function, which is sensitive to the matter 
density and expansion history of the universe, and to 
modifications of the laws of gravity. Measurements of 
the evolution of the galaxy cluster mass function across 
a broad range of redshifts can thus provide a powerful 
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tool for constraining the cosmological parameters (e.g., 
Vikhlinin et al. 2009a, 2009b) . Numerous galaxy cluster 
surveys will soon begin delivering a huge amount of data 
at optical. X-ray, and millimeter wavelengths, e.g. from 
Subaru/Hyper-Suprime-Cam, eROSITA, SPT and ACT. 
One of the main goals of these surveys is to measure the 
evolution of the galaxy cluster mass function, and thus 
to probe the expansion history of the universe. However, 
the mass of a galaxy cluster is not directly measurable. 
These surveys will therefore rely on "mass-like" observ- 
ables (e.g.. X-ray temperature - Evrard et al. 1996) and 
scaling relations between these observables and mass, to 
construct the all-important mass functions. Calibration 
of mass-observable scaling relations is therefore currently 
a high priority observational goal. 

Traditionally, observational studies of the mass- 
observable scaling relations have relied solely on X- 
ray observations, typically concentrating on the mass- 
temperature relation (e.g., Finoguenov et al., 2001; 
Sanderson et al., 2003; Ettori et al., 2004; Arnaud et 
al., 2005). X-ray-based mass measurements require hy- 
drostatic equilibrium (H.E.) and spherical symmetry to 
be assumed, and either measurement of the temperature 
profile, or an assumption of isothermality. Inclusion of 
X-ray temperature information in both axes of the mass- 
temperature relation may therefore induce intrinsic cor- 
relations into the measured relation. The validity of the 
underlying assumptions also warrants careful testing. 

Gravitational lensing offers cluster mass measurements 
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that are independent of X-ray observations, and do not 
rely on assuming H.E.. Joint lensing/X-ray studies (e.g., 
Okabe & Umetsu 2008; Kawaharada et al. 2010) are 
therefore a promising route for calibrating cluster mass- 
observable scaling relations. Indeed, early Icnsing/X-ray 
studies of cluster cores indicated that the scatter in clus- 
ter temperature may be as large as 40% at fixed mass, 
and that the scatter is dominated by disturbed, merging 
clusters, in which H.E. may not hold (Smith et al. 2005). 
Subsequent work has concentrated on using weak-lensing 
data to extend this pioneering work beyond cluster cores 
to overdcnsitics of 500 <J A ^ 2500 with respect to the 
critical density (Bardeau et al. 2007; Hoekstra 2007; Ped- 
ersen & Dahle 2007; Zhang et al. 2007, 2008). The main 
limiting factors in these weak-lensing/X-ray studies have 
been the limited statistical precision and heterogeneity 
of the available wcak-lcnsing data, and also the small 
samples observed to date. 

On the theoretical side, Kravtsov et al. (2006) 
proposed the so-called quasi-integrated pressure, 
Yx=Mgs,sXT as a "new robust low-scatter X-ray mass 
indicator" , or, a mass-like observable. This was mo- 
tivated by analysis of their hydrodynamic numerical 
simulations of clusters using an adaptive mesh refine- 
ment (AMR) code. They found that the temperature 
deviations from the M — T relation are anti-correlated 
with the gas mass deviations from the M — Mgas relation. 
This anti-correlation found in their simulations acts to 
suppress the scatter in the M — Ix relation, independent 
of the dynamical state of the clusters. This prediction 
has stimulated much observational effort within the 
X-ray community that has broadly supported the idea 
that Ix is the optimal X-ray mass proxy (e.g., Maughan 
2007; Arnaud ct al., 2007; Vikhlinin ct al. 2009a). 

However, Stanek et al.'s (2010) smoothed particle hy- 
drodynamic (SPH) Millennium Gas Simulations contra- 
dict Kravtsov et al.'s simulations. Stanek et al. predict 
that the temperature and gas-mass deviations are posi- 
tively correlated; this result appears to be independent of 
the range of gas physics (gravity only, cooling, preheat- 
ing) implemented in the simulations. Juett et al. (2010) 
have also recently suggested that previous X-ray-only 
studies may have underestimated the scatter in mass- 
observable scaling relations by a factor of '^ 2— 3. In sum- 
mary, a joint Icnsing/X-ray observational investigation of 
the relationships between mass and gas mass, tempera- 
ture, and Yx, is urgently needed. Such joint studies also 
lend themselves well to the task of observationally testing 
various corrections that have been derived from numeri- 
cal simulations to account for deviations from H.E.. For 
example, numerous authors have pointed out that H.E. 
mass estimates may underestimate the cluster mass be- 
cause of non-thermal pressure support due to turbulence 
caused by bulk motion of the cluster gas (e.g., Evrard 
1990; Rasia ct al. 2006; Nagai et al. 2007; Piffaretti & 
Valdarnini 2008; Fang et al. 2009), and Vikhlinin et al. 
(2009a) applied a 17% upward correction to X-ray masses 
of disturbed clusters, based on the results of simulations. 

A key goal of the Local Cluster Substructure Survey 
(LoCuSS ) is to calibrate cluster mass-observable scaling 
relations for future cosmological experiments. LoCuSS 
is a multi- wavelength survey of galaxy clusters at 0.15 < 

^|http://www.sr.bham.ac.uk/locuss| 



z < 0.3 selected from the ROSAT All-sky Survey cata- 
logs (Ebeling et al. 1998, 2000; Bohringer et al. 2004). 
To date we have published the first lensing/Sunyaev- 
Zeldovich effect comparison (Marrone et al. 2009), be- 
gun our Icnsing/X-ray scaling relation work with a pilot 
study (Zhang ct al. 2008), and compared lensing-based 
masses with H.E. masses on both small (Richard ct al., 
2010) and large (Zhang ct al. 2010) scales. This article 
is a continuation of our pilot study (Zhang et al. 2008), 
in which we combined weak-lensing mass measurements 
from the Canada-France-Hawaii Telescope (Bardeau et 
al. 2005, 2007) and from the Nordic Optical Telescope 
and UH SSin (Dahle 2006) with XMM-Newton observa- 
tions to calibrate the mass-observable scaling relations. 
As alluded to above, Zhang et al.'s results were lim- 
ited by the quality of the weak-lensing mass measure- 
ments, because the underlying data were heterogeneous 
in observing facilities, fields of view, and filters used. 
In this article we address these issues by using our own 
weak-lensing mass measurements based on uniform anal- 
ysis of our Subaru/Suprimc-cam observations (Okabe & 
Umetsu 2008; Okabe et al., 2010). Nevertheless, our 
Svth&in/ XMM-Newton sample remains small, at just 12 
clusters. As we discuss throughout this article, sample 
size therefore remains an issue, and we will address this 
in a future article. 

The outline of this paper is as follows. In Sec. [5] we 
briefiy describe the weak lensing and X-ray analysis, and 
measure the dynamical state of each cluster using XMM- 
Newton data. We present the main results on the mass- 
observable scaling relations in Sec. [3J discuss the results 
in Sec. SI and summarize our work in Sec. [S] Throughout 
this paper, we assume 0.mfi = 0.3, ^a = 0.7, and h ~ 
iJo/100 kms-^Mpc-i = 0.7. 

2. SAMPLE AND DATA ANALYSIS 

2.1. Sample 

For the purpose of this paper, we compiled a sample 
of 12 clusters -A 68, A 115, A 209, A 267, A 383, A 1835, 
A 1914, Z7160, A 2261, RX J2129.6-H0005, A 2390, and 
A 2631 - that represents the overlap between the sam- 
ples for which Subaru/Suprime-Cam and XMM-Newton 
data are available, and that we have previously pub- 
hshed (Zhang et al. 2008; Okabe & Umetsu 2008; Ok- 
abe et al. 2010). The sample does not suffer, by de- 
sign, any strong biases to extreme merging or extreme 
cool core clusters, and therefore can be regarded, qual- 
itatively, as representative of massive. X-ray luminous 
clusters. However, given the small sample size, we re- 
frain from attempting to quantify how these 12 might 
be biased with respect to the underlying cluster popu- 
lation in this article. Instead, this article presents some 
early results from our Suharu/ XMM-Newton program, 
that benefit from the use of our Subaru data, as opposed 
to the CFH12k/UH8k/NOT data that we used in Zhang 
et al. (2008). We defer detailed discussion of sample def- 
inition and possible biases to future articles in this series 
that will address larger, more complete samples. 

2.2. Weak-lensing mass measurements 

The details of our weak-lensing analysis are described 
in by Okabe & Umetsu (2008), and Okabe et al. (2010); 
here we provide a brief outline of some important aspects 
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of our methods. 

We selected background galaxies based on their loca- 
tion in the color-magnitude plane - typically {V — i')/i' - 
bluer or redder than cluster red-sequence by a minimum 
color-offset (Umetsu & Broadhurst 2008; Umetsu et al. 
2009; Okabe et al. 2010). As demonstrated by Okabe 
et al. (2010), contamination of the background galaxy 
catalogs by faint (unlensed) cluster members dilutes the 
wcak-lcnsing signal. This effect is more pronounced at 
smaller clustercentric radii because the number density of 
cluster galaxies rises towards the cluster centers. In the 
absence of our color-selection techniques, weak-lensing 
-^^500 E^nd A'/25oo measurements can be biased low by 

- 20%-50%. 

We used the COSMOS photometric redshift catalog 
(Ilbert et al. 2009) to estimate the redshift of the back- 
ground galaxies. Specifically, we calculated the average 
lensing weight, (I?ls/-Dos) = j^^ dzdPwh/dzDLs/Dos 
(see also Equation (10) in Okabe et al. 2010), of each 
background galaxy catalog by selecting galaxies identi- 
cal to both our catalogs, and the COSMOS catalog. Z?os 
and Dls are the angular diameter distances between the 
observer and source (background galaxy) and lens and 
source respectively. 

In cosmology the three-dimensional spherical mass. 
Ma, enclosed within a sphere of radius ta for a given 
overdensity A is most relevant for the cluster mass func- 
tion, where ta is chosen such that the average density 
within the sphere is equal to the critical mass density at 
the cluster redshift, pcr, times the overdensity A. We 
estimated Ma for each cluster by fitting the measured 
radial profile of lensing distortion signals to the NFW 
model prediction parameterized by the mass Ma and ca , 
where the NFW mass profile (Navarro, Frenk & White 
1996, 1997) is given as p ex r~^(l + cat/ta)^^ with ca 
being the concentration parameter. 

Describing cluster-scale dark matter halos as spherical 
objects may cause systematic errors in individual mass 
measurements because clusters are predicted to be triax- 
ial in the coUisionlcss CDM model (Jing & Suto 2002). 
For example, if the major axis of a triaxial halo is aligned 
with or perpendicular to the line of sight, a spherical 
model would overestimate or underestimate the mass, re- 
spectively, and also cause systematic errors in the mea- 
surement of the concentration parameter (Oguri et al. 
2005; Gavazzi 2005; Corless et al. 2009). However, if 
the distribution of cluster orientations is random, then 
adopting spherical mass models should not introduce a 
significant bias into the properties of the sample. We 
therefore check that this is the case for our sample by 
comparing the spherical mass measurements from Okabe 
et al. (2010) that we use here with triaxial mass measure- 
ments of the same clusters using the same background 
galaxy catalogs from Oguri et al. (2010). On average the 

spherical {M^^ ) and triaxial {M^" ) masses agree well 

- (Mi'"VA^A™^) = 0-98 ± 0.15,0.90 ± 0.17 and 0.83 ± 
0.21, for A = 500, 1000 and 2500 - confirming the expec- 
tation of negligible bias. 

2.3. X-ray observables 

The observations and data reduction are described in 
detail by Zhang et al. (2007, 2008). In brief, the three 
mass proxies considered in this article are calculated as 



follows. The global temperature is a volume average 
of the spectrally measured, radial temperature profile 
limited to the radial range of (0.2 — 0.5)r5oo. The gas 
mass Mgas('') was obtained for each cluster by integrat- 
ing a double-/? model of the electron density that was 
fitted to the X-ray surface brightness profile. The quasi- 
integrated pressure is the product of the gas mass and the 
global temperature: Yx{r) = Mgas(r) xTq. 2-0.57-500 ■ ^^te 
that Afgas(7'); To.2-o.5r5oo ^-^d Yx{r) have been calculated 
using radii obtained from the weak-lensing analysis, and 
not using radii calculated from the X-ray analysis as in 
Zhang et al. (2008). This definition of radii introduces 
a subtle correlation with weak lensing mass - we will 
explore this when estimating the intrinsic scatter in the 
mass-observable scaling relations in Sec. 13.61 and the Ap- 
pendix. Finally, we adopted a self-consistent definition 
of the cluster centers based on the weak lensing analysis. 
This caused us to change the centers of just two clusters 
- A 1914 and A 2631 - from those used by Zhang et al. 
(2010). 

2.4. X-ray morphology and dynamical state 

Previous joint lensing/X-ray studies have identified the 
dynamical state of clusters as a significant source of scat- 
ter in mass-observable scaling relations (Smith et al. 
2005; Pedersen & Dahle 2007; Zhang et al. 2008, 2009). 
In this section we therefore classify the clusters as either 
"disturbed" or "undisturbed", based on a new method 
patterned on those developed for the morphological clas- 
sification of galaxies (e.g., Conselice 2003). 

We calculate the asymmetry [A] and fluctuation (F) 
of the X-ray surface brightness distribution in the 0.7 — 
2 keV band. Asymmetry is defined as ^ = iJ^a l-^ij ~ 
Rij\)/J2ii ^'J' ^^^ normalized sum of the absolute value 
of the flux residuals where Jy is a matrix element of the 
combined M0SH-M0S2 XMM-Newton frame in the 0.7- 
2.0 kcV band, flat flelded, point source subtracted and 
refilled assuming a Poisson distribution, and Rij are the 
matrix elements obtained by rotating the above frame by 
180°. The pixel size of both frames is 4" x 4". The fluc- 
tuation, F, measures deviations from a smooth flux dis- 
tribution and is defined as F = {J^ij hj - Bij)/ J^ij hj, 
where Bij is an element in a frame smoothed on 2 arcmin 
scales, which corresponds to a physical scale of 400 kpc at 
z — 0.2. Such smoothing also suppresses the effect of the 
complex shape of the XMM-Newton point spread func- 
tion (Ghizzardi 2001). We estimate the statistical errors 
of A and F assuming Poisson noise computed within a 
radius of r^oQ, excluding CCD gaps and bad pixels. We 
also estimate the systematic error of A caused by uncer- 
tainties in the cluster centers by recalculating A, each 
time moving the cluster centers onto one of the neigh- 
boring pixels within the r < 4" circle from the nominal 
cluster center. 

The clusters span the range A ^ 0.07 — 0.15 and 
F ~ — 0.14 (Fig. [T|). Dynamically disturbed clusters 
generally have an asymmetric X-ray morphology, with an 
offset between optical and X-ray centers, and arc there- 
fore expected to have larger A and F than undisturbed 
clusters. To separate the clusters into two subsamples 
that represent relatively disturbed and relatively undis- 
turbed systems, we subdivided the A — F plane into four 
quadrants: (1) A < 1.1 and F < 0.05 - RX J2129, A209, 
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A 383, A 1835, and A 2390, (2) A > 1.1 and F < 0.05 
- A 2261 and A 1914, (3) A < 1.1 and F > 0.05 - 
A 68, A 2631, A 267, and Z7160, and (4) A > 1.1 and 
F > 0.05 - A 115. We classify the five clusters m quad- 
rant (1) - low A and low F - as undisturbed clusters, 
and the remaining seven as disturbed clusters. It is im- 
mediately obvious that this classification matches other 
possible classification schemes well. For example, four of 
the five undisturbed clusters host a cool core (e.g. Smith 
ct al. 2003; Allen et al. 2001; Peterson et al. 2002), and 
the disturbed clusters have been discussed extensively as 
merging/cold- front clusters (e.g., Okabe & Umetsu 2008; 
Mazzotta & Giacintucci 2008; Gutierrez & Krawczynski 
2005), in which complicated temperature/entropy distri- 
butions or large offsets between lensing/optical and X- 
ray centroids exits (e.g., Finoguenov et al. 2005; Smith 
et al. 2005; Sanderson et al. 2009a). In summary, all of 
the clusters identified as disturbed in the A — F plane arc 
independently confirmed as disturbed by other methods 
in the literature. However we stress again the relative na- 
ture of the disturbed/undisturbed classification, and ac- 
knowledge that the disturbed clusters in particular likely 
comprise clusters in a wide variety of stages in their dy- 
namical evolution. We will return to this issue later when 
we assess the impact of a single extreme merging cluster 
on our attempts to calibrate the mass-observable scaling 
relations. 

3. RESULTS 

In this section, we present the main empirical results 
of the slope, normalization, and intrinsic scatter in the 
mass-observable scaling relations and how these depend 
on the dynamical state of the clusters. We also discuss 
the correlation between gas mass and temperature devi- 
ations. 

3.1. Scaling relations and fitting methods 

If gravitational heating is the dominant mechanism re- 
sponsible for the X-ray properties of galaxy clusters, the 
following scaling relations are expected to hold: 



ME{z) ex AfgasF(z)/i3/2^ 
ME{z) oc T^^^h^\ 



(1) 
(2) 
(3) 



where M, Mgas, and T are the total mass, gas mass, and 
temperature of a cluster, respectively, and Yx = Mgas xT 
is the quasi-integrated pressure. These relations, specif- 
ically the exponents of M, Mgas, and T, are usually re- 
ferred to as self-similar, following Kaiser (1986). Note 

that the term E{z) = H{z)/Hq = [rj„,o(l + 2)^ + ^a]^/^ 
accounts for the rcdshift evolution of the clusters in a flat 
universe. 

In the following subsections we therefore fit the func- 
tional form M^ = MqX2 to the data, where M^ = 
M E{z), Mq is the normalization, Xz is the X-ray observ- 
able (i.e., Yx, T, or Mgas) multiplied by E{z) or not as ap- 
propriate based on Equations 1-3, and 7 is the logarith- 
mic slope. These fits are done at three overdensities with 
respect to the critical density: A = 2500, 1000 and 500. 
The scaling relation slope and normalization measure- 
ments are based on orthogonal regression performed us- 
ing the Orthogonal Distance Regression package (ODR- 



PACK, e.g. Boggs et al. 1987) taking into account the 
measurement errors. In general we ignore the subtle 
correlations introduced by measuring the X-ray observ- 
ables within radii defined by the weak-lensing analysis, 
although we do take them into account in Sec l3.3l when 
we measure the intrinsic scatter. To check for consistency 
with other work, we have also refitted the relations using 
the bisector modification of the BCES method (Akritas 
& Bcrshady 1996). The difference on best-fit scaling re- 
lation parameters between the two fitting methods is a 
small fraction of statistical uncertainties. For example, 
the difference on best-fit slopes and normalizations be- 
tween the two methods is typically ^ 30% and ^ 6% of 
the statistical error respectively. We also did the boot- 
strap resampling to estimate the sample variance on the 
slope parameter, and found that it is <^ 20% of the sta- 
tistical errors. 

3.2. Slope and normalization 

We first fit the scaling relations to the full sample of 12 
clusters with both slope 7 and normalization Mq as free 
parameters. At A = 500 the best-fit slopes of all three 
relations agree well with the self-similar model (Tabled)). 
At higher overdensities, the agreement deteriorates for all 
three relations, indeed the slopes of the M 



2500 



M, 



M2500 ~ Yx relations are discrepant from self-similar at 
~ 2 — 3(7 at A = 2500. This fiattening in the scaling 
relations at higher A can also be seen graphically in Fig- 
ures H [3] and El in which we show the Ma-Ix, Ma-Tx 
and Ma -Mgas relations respectively. 

To constrain the normalization parameter Mq, we fix 
the slope parameters to the self-similar values and repeat 
the fits. The measured normalizations are all consistent 
with those obtained by Zhang et al. (2008) using the 
same XMM-Newton data and independent weak-lensing 
data. The superior quality and uniformity of our Sub- 
aru data shows differences between the normalizations 
for disturbed and undisturbed clusters. These differences 
are most pronounced at A = 500 - see Table [2] ~ specifi- 
cally, at fixed Yx undisturbed clusters are measured to be 
^ 22% more massive than disturbed clusters at ~ 1.5cr 
significance. Similarly, at fixed T undisturbed clusters 
are measured to be ^ 43% more massive than disturbed 
clusters at ^ l.Sa significance. We confirm that our re- 
sults are insensitive to whether or not the slopes are fixed 
to the self-similar value. 

3.3. Scatter 

We also measured the intrinsic scatter, (Tinj\/, for 
the logarithm of the F-axis, ME{z), for each mass- 
observable scaling relation using the Bayesian method 
described in the Appendix. Here we take into account 
the correlations caused by measuring X-ray observables 
within radii defined by the lensing analysis - see the 
Appendix. We also confirmed that the best-fit slopes 
and normalizations obtained using the ODR methods 
discussed above are consistent within errors with those 
obtained using the more sophisticated Bayesian method 
considered here. The intrinsic scatter in all three rela- 
tions is well described by a lognormal distribution. 

The M~T relation exhibits the largest intrinsic scatter 
(^ 0.23 — 0.33; Table[3]) among the three mass-observable 
relations. We also observe an increase in the intrinsic 
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scatter with increasing radius (i.e. decreasing the interior 
overdensity A ) . The same trend is found in undisturbed 
clusters, while the opposite trend is found in disturbed 
clusters. However, this trend is not a physical feature 
of the intracluster gas affected by gravitational heating, 
because we used a fixed global temperature measurement 
in the radial range of 0.2— O.Srsoo for all the overdensities. 

The M — Mgas relation is the tightest of the three, 
with an intrinsic scatter in mass of (J\nM ~ 0.12 — 0.16 at 
A = 500 and 1000. At A = 2500 the scatter is roughly 
double that at lower overdensity (Table [3|), which may 
be due to different core properties of individual clusters. 
For example, cool-core clusters have denser, cuspier cores 
than non-cool core clusters (e.g., Croston et al. 2008; 
Sanderson et al. 2009b; McCarthy et al. 2008). Such 
differences between cluster cores have a much smaller ef- 
fect on measurements at larger radii because the core 
regions make a small contribution to the total gas mass 
measured out to A = 1000 and 500. However, note that 
the intrinsic scatter is not well constrained for M — Mgas 
because the scatter is dominated by statistical errors. A 
larger sample is clearly needed to improve the constraints 
on the intrinsic scatter in M — Mgas, however, it is im- 
portant to note that this is the only relation that appears 
to have ■^ 10% intrinsic scatter. 

The observed intrinsic scatter in the M — Yx relation is 
intermediate between that of the M — Mgas and M — T 
relations, at CTinj\f ~ 0.20 — 0.25 (Table [3l). This is a 
factor of > 2 greater than that originally predicted by 
Kravtsov et al. (2006) based on their AMR simulations. 

3.4. The impact of an outlier 

In this section we highlight the impact of one cluster, 
A 1914, on our results. This cluster has previously been 
identified as a merging cluster with a complex X-ray mor- 
phology, radio halo, and weak-lensing-based dark matter 
distribution (Buote & Tsai 1996; Bacchi et al. 2003; Gov- 
oni et al. 2004; Okabe & Umetsu 2008). We have also 
identified it as having the most extreme X-ray/lensing 
mass discrepancy among the 12 clusters considered here 
(Zhang et al. 2010). 

To assess the impact of such clusters on the measured 
reliability of X-ray observables as mass proxies we re- 
peated the calculations of normalization and scatter dis- 
cussed in Sec. [3J]&p3] excluding A 1914 (Tabled. At 
A = 500 the normalization of the M — Yx relations for 
disturbed and undisturbed clusters are different at just 
^ 1.2cr significance when A 1914 is excluded from the dis- 
turbed sample, in contrast to the ~ 1.8ct difference based 
on the full sample of 12 clusters. We also find that ex- 
cluding A 1914 reduces the intrinsic scatter on all of the 
scaling relations. In particular, the intrinsic scatter on 
M - Yx is reduced by - 25% from ainM ~ 0.20 to 0.15. 
Jack-knife tests on samples of 11 clusters (i.e. removing 
each cluster in turn) also confirm that A 1914 is indeed 
the most significant outlier among our sample. 

These results indicate that outliers in the cluster pop- 
ulation require careful treatment in the construction and 
application of mass-observable scaling relations. In sum- 
mary, reliable cluster selection functions are required to 
gain robust constraints. This will be especially true for 
future high redshift surveys because the fraction of merg- 
ing clusters is expected to increase with look-back time 



(Vikhlininet al. 2009a). 

3.5. Covariance of deviations 

We investigate the covariance of deviations from the 
best-fit M — Mgas and M — T relations following re- 
cent numerical simulation studies (e.g. Kravtsov et al. 
2006; Stanek et al. 2010). For a given mean scaling rela- 
tion Y = f{X), the deviations of each cluster from the 
mean relation are quantified a.s 5Y = [Y — f{X)] and 
SX = [X — f^^{Y)]. We use the mean normalizations 
for a full sample of 12 clusters, however we found that 
the following results do not change significantly when 
the best-fit normalizations of undisturbed and disturbed 
clusters are used instead. 

The temperature and Afgas deviations, Sr/T{M^) 
and i5Mgas/Mgas(MA), appear to be positively corre- 
lated (Fig. [B]). We test this quantitatively using Spear- 
man's rank correlation coefficient test, obtaining r^ ~ 
0.531 ± 0.009. The probability of obtaining a value of 
Vs greater than or equal to the measured value is low: 
V = 0.075 ±0.006. This test therefore indicates that the 
positive correlation is significant. However, the appar- 
ent positive correlation between the temperature and gas 
mass deviations does not show the correlation between 
intrinsic scatter, but between total scatter which is a 
convolution of measurement errors and intrinsic scatter, 
because we here did not take into account for measure- 
ment uncertainties. When dealing with observational 
constraints on scaling relations, it therefore essential to 
include both the covariance of intrinsic scatter, and the 
measurement errors with which the scatter is convolved 
in robust calculations. 

3.6. Covariance of intrinsic scatter 

We simultaneously fit M — T and M — Mgas relations 
and measure the covariance of the intrinsic scatter us- 
ing a multi-dimensional fitting method described in the 
Appendix. This method considers not only the matrix 
of the observational errors for individual clusters, Sobs,i, 
but also the covariance matrix of the intrinsic scatter, 
Sint. The covariance of the intrinsic scatter is given by 



(Jtg 
2 



<^tq cr. 



(4) 



where a^ and a^ are variances for the logarithm of tem- 
perature and gas mass, respectively, and atg = rutOg is a 
covariance with a coefficient r. Here, we do not need to 
take into account for intrinsic scatter on mass, because 
the gas properties, under a cluster mass given by the cos- 
mology, only physically have intrinsic scatter due to the 
gas evolution. When we estimate cluster masses from 
X-ray observables via scaling relations, there is intrinsic 
scatter on mass due to the propagation from gas intrin- 
sic scatter. As shown in the Appendix, the observational 
error matrix for individual clusters is given by 



Ki 



9 n. 
d In M, 



d In M 
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where e™ , et and Cg are observational errors for the log- 
arithm of the mass, temperature and gas mass, respec- 
tively. The coefficients of e™ are the slopes of the mass 
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observable relations T oc M^/'' and Mgas oc M, and a 
term to account for propagation of the error on ta (de- 
rived from the lensing analysis) to the error on the X-ray 
observables, 91nMgas/91nM (see the Appendix). Since 
the correlation between the observational errors of tem- 
perature and gas mass is negligible, we set them to zero. 
We measure the intrinsic scatter for the M — Mgas and 
M — T relations {(Jg,<Jt and r) with their slopes set to 
the self-similar prediction. The likelihood (eq. IA2[) is 
given in the Appendix. We explore parameter space us- 
ing the Metropolis Hastings algorithm, as described in 
the Appendix, restricting ranges of values explored to 
< cTt < 1, < o-g < 1 and \r\ < 1. 

We first impose a flat prior on at,<7g, and r within 
the limits referred to above, and obtain measurements of 
intrinsic scatter that are in good agreement with those 
obtained by fitting the M — Aigas and M — T relations in- 
dependently: ag = 0.132j:o:o86 and at = 0.213l:o:o8o (Ta- 
bleU}. However, the coefficient, r, is not well determined 
due to the tail of the posterior probability distribution 
extending to negative values. Nevertheless, we derive a 
68.3% confidence lower limit of r ;> 0.185. We also mea- 
sure the mode of the marginalized posterior probability 
distribution to be r = 0.575. 

We then repeat the fit, this time using Gaussian pri- 
ors centered on the best-fit measurements of at and ag. 
As shown in the second column in Table |4l all resulting 
parameters are consistent with those in the flat prior. 
This time the mode of the posterior probability distri- 
bution is r = 0.570., and we find that the lower limit 
on r is again positive - i.e. it is positive independent of 
the prior. The positive coefficient leads to a large in- 
trinsic scatter of Yx 



ol + a? 



MgasT because the last term in 
2ratag becomes positive. The posi- 
tive coefficient indicates that deviations in gas mass and 
temperature are partially correlated. 

However, we stress that the small size of our observed 
sample severely limits the statistical power of our results 
- specifically, we cannot rule out the possibility that the 
gas mass and temperature deviations are two random 
variables that are correlated by accident. To quantify 
this, we calculated the probability, "Pdr'l > \r\), that the 
correlation coefficient of the two random variables in a 
sample of 12 drawings, r' , is higher than the observed 
value (Pugh & Winslow 1966). Since the probability is 
non-negligible (T'dr'l > |r|) < 0.565 and 0.597; Table g]), 
we cannot yet completely exclude this possibility. There- 
fore, we need to increase the sample size before we can 
make definitive statements on the correlation of intrinsic 
scatter. 

4. DISCUSSION 

4.1. Comparison with previous observations 

First we compare our results with those of our pi- 
lot study (Zhang et al. 2008). The main difference be- 
tween Zhang et al.'s analysis and that presented here is 
that Zhang et al.'s weak-lensing mass measurements were 
drawn from the literature, and thus suffered heteroge- 
neous image quality, observed depth, and systematic un- 
certainties relating to background galaxy selection and 
faint galaxy shape measurement. Despite these differ- 
ences, the overall normalization of our mass-observable 
scaling relations agree within the uncertainties with those 



of Zhang et al. (2008). However despite our sample be- 
ing roughly a factor of two smaller than that of Zhang et 
al., we detect structural segregation in the M — Yx and 
M—T relations at ~ 2a significance. Our ability to make 
this detection is likely due to the factor of ^ 2 smaller 
statistical errors on weak-lcnsing mass measurements of 
individual clusters, thanks to the superb quality of our 
Subaru data. 

This is the first time that structural segregation has 
been found in the M — Yx relation, however it has been 
detected in the M — T relation at a similar amplitude, 
and level of significance in previous joint lensing/X-ray 
studies (Smith et al. 2005; Pedersen & Dahle 2007). The- 
oretical studies (e.g. Randall et al. 2002) suggest that 
this segregation may be caused by cluster-cluster mergers 
boosting the temperature of disturbed merging clusters. 
However orientation effects may also contribute to sys- 
tematic errors in weak-lcnsing mass measurements that 
exaggerate the segregation (Meneghetti et al. 2010). We 
therefore defer physical interpretation of the observed 
segregation to a future careful investigation of the degen- 
eracy between cluster orientation and residuals on mass- 
observable scaling relations. 

We also compare the normalizations of our M — Yx and 
M — T relations at A = 500 with the same from Vikhlinin 
et al.'s (2009a) X-ray-only study of "relaxed" clusters 
with Chandra. Note that, as discussed in Sec. l3.1[ our re- 
sults are insensitive to whether we use orthogonal regres- 
sion or bisector fitting techniques; Vikhlinin et al. used 
the bisector method. The normalizations of our M — Yx 
and M — T relations for undisturbed clusters agree within 
the uncertainties with Vikhlinin et al.'s (2009a) relaxed 
clusters (Figs. [2] & |3]). This suggests that orientation 
effects may not be a major influence on the normaliza- 
tion of our undisturbed cluster scaling relations. On the 
other hand, it is thus clear that the normalizations of our 
disturbed cluster scaling relations differ from Vikhlinin 
et al.'s results at ~ 2ct significance. It will therefore 
be important in the future to study the mass-observable 
scaling relations with large samples comprising clusters 
that span a wide range of dynamical states. This need 
will become more acute at high redshift because the frac- 
tion of disturbed clusters likely increases toward higher 
redshifts. 

4.2. Comparison with simulations 

The normalizations of our M — A/gas and M — Yx 
relations arc lower in mass at fixed X-ray observable 
(Fig. [2] & [3]) than the predictions from Nagai et al.'s 
(2007) simulations. This result is consistent with our 
pilot study (Zhang et al. 2008), in which we compared 
our joint lensing/X-ray observational results with a wider 
range of simulations including those of Borgani et al. 
(2004) and showed that the normalizations disagree with 
the simulations at > 2a for M — Yx and > Icr for Af — T 
and M — A'/gas relations, respectively. 

The structural segregation found in our M — T and 
M — Yx relations - i.e. undisturbed clusters are ~ 40% 
and ^ 20% more massive than disturbed clusters at fixed 
temperature and fixed Yx respectively at ^^ 2a signifi- 
cance - disagrees with Kravtsov et al.'s (2006) simula- 
tions, upon which they based their proposal of Yx as 
a low-scatter mass proxy. In summary, they found no 
difference in the normalization of M — Yx between clus- 
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ters that they classified as "relaxed" and "unrelaxed" 
(roughly equivalent to our undisturbed/disturbed classi- 
fication), and that their relaxed clusters are less massive 
at fixed temperature than unrelaxed clusters - i.e. oppo- 
site to our observational result. 

The intrinsic scatter in the M — T and M — Mgas rela- 
tions at A := 500 (~ 20-30% and ~ 10% respectively) is 
comparable to Kravtsov et al. (2006). However the anti- 
correlation between temperature and gas mass deviations 
from the mean M — T and M — Mgas relations predicted 
by Kravtsov et al., and that has motivated much atten- 
tion on Yx as a low scatter mass proxy, is not supported 
by our data, although we need to investigate this for a 
larger sample. In contrast with Kravtsov et al. (2006), 
we find that the lower limit of the positive coefficient 
for intrinsic scatter is in better agreement with Stanek 
et al.'s SPH Millennium Gas Simulations showing that 
the coefficients between the spectroscopic-like tempera- 
ture and gas fraction for z = and A ~ 200 are positive 
irrespective of the process of preheating and cooling. 

Correlated gas mass and temperature deviations imply 
a possibility that adiabatic compression/expansion of the 
intracluster gas is important in cluster evolution. These 
adiabatic fluctuations propagate much faster than cool- 
ing losses, and thus help to explain why Stanek et al.'s 
simple hydrodynamical simulation matches the observa- 
tions well. Umetsu et al. (2010) have also suggested that, 
in the adiabatic expansion phase of a post-merger, both 
temperature and (encompassed) gas mass decrease com- 
pared to those before the merger. Another possibility 
can be due to the departure from the spherical symme- 
try. As long as the intracluster gas is in H.E., even in 
the elongated gravitational potential well, the gas mass 
distribution is supposed to be more round than the dis- 
tribution of dark matter. In our modeling, we assume 
spherical symmetry of the gas and mass proflles, in which 
deviations from spherical geometry are transferred into 
systematic measurement and may in part lead to a pos- 
itive correlation between the gas mass and mass. Fur- 
thermore, many possible deviations in mass-observable 
relations could be relevant to the interpretation of these 
results. For example, star-formation efficiency affects the 
total gas mass and the epoch for cluster formation af- 
fects the cluster temperature. Some of these effects imply 
other correlations that can be explored in the future with 
a larger sample, e.g. correlations between temperature 
deviations and dark matter profile shapes, or between 
deviations in the gas mass and total stellar fraction. 

4.3. Principal component analysis of observational 

data 

Finally, we propose a method, based on principal 
component analysis, for constructing and calibrate a 
low-scatter mass proxy using solely observational data. 
Thanks to multi-dimensional fitting, since we obtained 
both the intrinsic covariance and the normalization, we 
here do not need to take into account for measurement er- 
rors. By solving the eigenvalue equation (S — cr^J)y ~ 0, 
where I is the identity matrix, we can obtain its eigen- 
values CTj. and eigenvectors y± = \nY±, as follows: 



P±' 



<7l±{{a 



^2\2 



Ar'a^alf' 



Y±=Mg^TP^, 



(5) 
(6) 



2rat(7„ 



1/2 



(7) 



where Y- is the mass proxy with the smaller scatter (T_ . 
Note that if the coefficient r is negative, as predicted by 
Kravtsov et al. (2006), then the temperature exponent, 
P-, is always positive. In this framework, the scaling 
relation can be written as: 

ME{z) OC (y±£;(z))'/('+'P±) ;^(9-4p±)/(3+2p±)/2^ (g) 

and the intrinsic scatter on mass when one estimate from 
X-ray observables via new scaling relation is given by: 

3 



ClnA/ — 



|3 + 2p± 



■a±. 



(9) 



Basically, a combination of highly correlated/anti- 
correlated observables gives a new mass proxy with 
smaller scatter. Unfortunately, our sample is too small 
to constrain the covariance of the intrinsic scatter well, 
and to perform this principal component analysis. This 
exercise awaits enlargement of our sample. In princi- 
ple, a combination of the principal component analysis 
and the method for measuring the covariance of intrinsic 
scatter (see the Appendix) could be applied to multi- 
dimensional (i.e. weak lensing. X-ray, SZ and optical ob- 
servables) data sets both within LoCuSS and in large 
forthcoming surveys. 

5. SUMMARY 

We have presented a joint weak- lensing and X-ray anal- 
ysis of 12 clusters, based on Subaru and XMM-Newton 
observations as part of the LoCuSS Survery. The main 
goal is to calibrate the scaling relations between clus- 
ter mass obtained from weak lensing observations (M), 
and X-ray observables, specifically the gas temperature, 
gas mass, and quasi-integrated pressure (T, Mgas, and 
Yx = Mgas X T). An accurate understanding of these 
relations will be essential to the success of future at- 
tempts to constrain dark energy with clusters via growth 
of structure experiments. Our main results are summa- 
rized below. 

• The dynamical state of clusters can be diagnosed 
empirically via morphological classification of clus- 
ters based on asymmetry (A) and ffuctuation (F) 
parameters derived from imaging data from XMM- 
Newton. Undisturbed clusters are identified as 
those with relatively symmetric and smooth X- 
ray morphology - A < 1.1 and F < 0.05; dis- 
turbed clusters satisfy either or both of these cri- 
teria. Five clusters are classified as undisturbed 
and seven as disturbed. This classification matches 
those based on alternative measures such as the 
presence/ absence of cool cores, cold fronts, and 
substructures in lensing mass maps. 

• We detected structural segregation in the Al — T 
and M — Yx relations at A ~ 500, in the sense that 
undisturbed clusters are ^ 40% and ~ 20% more 
massive than disturbed clusters at fixed T and Yx 
respectively, at ~ 2a significance. Segregation in 
the M — T plane is qualitatively in agreement with 
some of the previous observational results (Smith et 
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al. 2005; Pedersen & Dahlc 2007); as far as we know 
these are the first joint lensing/X-ray results on 
the M — Ix relation. These results both contradict 
Kravtsov et al.'s predictions upon which they based 
their proposal that Yx may be a useful low-scatter 
mass proxy. 

• The intrinsic scatter in the observed M—T, M—Yx 



and M — Afgas relations is measured to be cti 



nM 



gas 



0.3, ~ 0.2, and - 0.1 respectively at A = 500. M^ 
therefore appears to be the most promising mass 
proxy of these three observables, especially because 
the scatter in M — Mgas appears to be independent 
of the cluster dynamical state. 

• The best-fit mass-observable scaling relations are 
sensitive to the inclusion/exclusion of one cluster in 
our sample, namely A 1914, a well known merging 
cluster. If this cluster is excluded from our analysis 
then the scatter is greatly reduced - most notably, 
the intrinsic scatter in the M — Yx relation is re- 
duced by ~ 25%. We conclude that a larger, more 
complete sample of clusters is required to reliably 
calibrate the scaling relations and to robustly mea- 
sure how the most extreme merging clusters influ- 
ence the relations. This will be particularly valu- 
able as scaling relation studies proceed to higher 
rcdshifts at which merging clusters are expected to 
become more prevalent. 

• Temperature deviations from the best-fit M — T re- 
lation and gas mass deviations from the M — Mgas 
relation are positively correlated. The coefficient 
between the gas mass and temperature deviations 
is positive, independent of our analysis methods, 
and is found to be r > 0.185. This result, in 
particular the lower limit on r agrees well with 
predictions based on Millennium Gas Simulations 
(Stanek et al. 2010), and disagrees with predic- 
tions based on Kravtsov et al.'s (2006) simulations. 
However, we caution that the chance probability, 
Pdr'l > |r|), that the correlation coefficient of the 
two random variables in a sample of 12 drawings, 
r', is higher than the observed one, is not small. We 
therefore cannot exclude the possibility that the 
measured correlation between gas mass and tem- 
perature residuals is an accident. A larger sample 
of clusters is needed to achieve definitive results. 

• Finally, we outlined a new method for constructing 
a robust low-scatter mass proxy, Y_ = MgasT^^, 
calibrated solely by observational data, based on a 
principal component analysis. This is a generaliza- 
tion of the quasi-integrated pressure Yx = MgasT, 
proposed by Kravtsov et al. (2006). In principle, 
a combination of the principal component anal- 



ysis and the method for measuring the covari- 
ance of intrinsic scatter could be applied to multi- 
dimensional data sets in order to construct robust 
new mass proxy. 



Our future program will concentrate on expanding the 
sample of clusters for which high quality Subaru and 
XMM-Newton data are available, in order to achieve 
definitive results on the issues raised in this article. 
Key issues will include improvement of the statistical 
uncertainties on the scaling relation fits, characterizing 
more fully the influence of extreme merging clusters on 
the scaling relations, and exploring the balance between 
physical and orientation effects in causing the observed 
structural segregation in the M—T and M — Yx relations. 
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APPENDIX 
MULTI-DIMENSIONAL FITTING OF DATA WITH COVARIANCE OF INTRINSIC SCATTER 

The measurement of a covariance of intrinsic scatter beyond observational errors is of prime importance in order to 
understand their intrinsic characteristics. Wc derive the multi-dimensional fitting with the covariance of intrinsic scat- 
ter, in the context of the Markov Chain Monte Carlo (MCMC) method with standard Metropolis-Hastings sampling, 
taking into account observational errors. 

Here, we suppose a data set (x = {xi,yii,y2i, ■ ■ ■ ,ypj}"=o) ^^ "■ sampling numbers and p + 1 variables. A linear 
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regression equation defined by Op + bpXi and an intrinsic covariance matrix Cint (p x p) in the y-coordinates are applied 
to modeling the relationship between x and t/p in the data set. The diagonal elements in the covariance of intrinsic 
scatter, Cint are intrinsic scatter, a?,, for the yp observables. The off-diagonal elements describe the intrinsic covariance 
between yp variables. 

A fitting parameter 6 is composed of the linear regression parameters, a, 6, and the elements of the intrinsic covariance 
matrix. From Bayesian statistics, the posterior probability of the parameter vector 6 is proportional to the conditional 
probability of x, p{x\6) and a prior probability function, Pprior(^), 



p{6\x) (Xp{x\e)pprior{0). 

The conditional probability of x, given parameters 6, is the likelihood described by 

1 



Pix\e) = n 



\ (27r)p/2|det(Ci)|i/2 



exp [-{vi - (a + bx)fCi ^{vi - (a + bx))/2] 



where the covariance matrix Ci = Ci^hs + Cint with the observational error covariance matrix, 
covariance. The loglikelihood of eq. (jA2p is given by 

-2L = Y, log(det(Ci)) + Y^iVi -{a + bx)fCi-\yi - [a + bx)). 



^dht 



(Al) 



(A2) 



and intrinsic 



(A3) 



The second term on the right-hand side of Equation (|A3|) is referred to as the y^ in some papers (e.g. Akritas & 
Bershady 1996; Tremaine et al. 2002; Pizagno et al. 2005; Novak et al. 2006; Weiner et al. 2006), and the intrinsic 
scatter is estimated by requiring the reduced x^ to be unity, in the framework of the x^ minimization fitting. It is, 
however, an inadequate fitting procedure because the first term, depending on the parameters, cannot be ignored (see 
also D'Agostini 2005). We therefore employ the likelihood function for a calculation of the covariance of intrinsic 
scatter. The diagonal elements in the observed error matrix are given by {hp — dyp^i/dxi)^e1 , -I- ef, , with the observed 



variances ei and e^, in x and 

J- Up 



"Vp^^ 



r,{x) variables. Here, dyi/dxi represents the error propagation between two variables. 



The off-diagonal elements, Cppi^i, for i sample are {hp — dyp^i/dxi){bpi ~ dyp'j/dxi)e^j + {e 



i). It describes the 



correlation between observational errors : the first term is an error correlation via the same x values, and the second 
one is the correlation between observational errors of yp and ypi . In this study, we measure the covariance of intrinsic 
scatter in the parameter plane of x = In(Af), y ~ (ln(T), In(Mgas)). 
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TABLE 1 
Slopes of the mass-observable relations for the full sample. 



Relation 
(1) 


7 using A/wL 
A : 500 1000 2500 

(2) (3) (4) 


Ma-Y^ 


0.67 ±0.14 0.59 ±0.11 0.46 ±0.11 


Ma-TZ,^o.s..,, 


1.49 ±0.58 1.49 ±0.46 1.26 ± 0.35 


Ma-MS^ 


0.98 ±0.15 0.86 ±0.14 0.68 ±0.15 



Note. — A single power law form of slope 7 is considered. The X-ray temperature is derived by the volume average of the deprojected 
radial temperature profile and using the radii within 0.2 — 0.5r5oo (T = Tb.2— o.Srsoo i ^^^ Zhang et al. 2008). Column (1): Scaling relations. 
Columns (2-4): Slopes of the mass-observable relations using weak Icnsing masses at the overdensities of A = 500, 1000 & 2500, respectively. 



TABLE 2 
Normalization (Mq) and morphological dependence of the mass-observable relations. 



Relation 
(1) 


All 12 clusters 

all undisturbed disturbed 
(2) (3) (4) 


w/o A1914 

all disturbed 

(5) (6) 


M500->X 

Miooo-Yx 
M2500-YX 


4-y^_0.33 0-*"-0.52 ^•*'-0.37 

4.5lj:«;f, 4.77l«:fo 4.28t«;f6 
4.02t»i? 4.06t«;f3 3.99t°;« 


5 16+"'^* 4 83+"'*^ 
4 70+0.28 . fi9+0.41 
*-'U_0.26 *-O^-0.37 
- r,r, + 0.30 . qo + 0.43 
4-20_o.2g 4.38_o.3g 


A^500-To.2-0.5r5oo 
A^1000^To.2-0.5r5oo 
-'^'^2500-^0.2-0. 51-500 


2 45+0-2''' 2 Q4+0-5U 2 05+"'2'* 
^•^0-0.24 ^-^^-0.42 ^-^0-0.22 
1 70+0-15 1 04+0-24 1 c-r|+0-18 
l-'U_oi4 l-y4_o.21 l-t'U_o 16 
(,07+0.08 1 f(o+0.10 PI 01+0.12 


2.34±i!:f2 2.2o±li:;i^ 

im+Jo\l 1.64±°;i? 


-MsoQ-Mgas 

-Wiooo-.^fgas 

-M2500-Mgas 


13.10+0;^^ l3.79+i;t^ 12.53+i;'^^ 

14.46tr8? 14.5211-Jl UAltlil 
17.33+li« 16.85+2-54 n.76+\-ll 


13.59±ii;'(^ 13.36±{;li« 

14.9lir84 15-36ll:S^ 
18.0lll;« 19.45lliJ 



Note. — The forms of M-Yx, M-Tx, and M-M^.^^ relations are given by MAE{zf/-' = Mo(Yx/3 X lOi4M0keV)5/3 x lO^^h^/'^MQ, 
MaE{z) = Mo{kBT/5 keV)3/2 ^ lOi^fe-lM©, and MaE(z) = Mo{MgasE{z))h^^^MQ, respectively. We fix the slopes to the self- 
similar values. The X-ray temperature is derived by the volume average of the radial temperature profile in the range of 0.2 — 0.5r5oo 
(T = 7b.2— o.5rroo ; see Zhang et al. 2008). Column (1): Scaling relations. Columns (2-4): Normalization of the fit to the relation using 
weak lensing masses. The results for the full sample, undisturbed and disturbed clusters are presented. Columns (5-6): Normalization of 
the mass-observable relations for the 11 clusters and 6 disturbed clusters, excluding A1914, respectively. 
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TABLE 3 

Intrinsic scatter in mass - observable relations 



Relation 

(1) 


All 12 clusters 

all undisturbed disturbed 

(2) (3) (4) 


w/o A1914 

all disturbed 
(5) (6) 


Miooo-yx 

A^2500->X 


0.203±;j;*^ < 0.283 0.2161^66 
0.203+[j-!^^| < 0.283 0.241+5;°^^ 
0.245t0;»g < 0.353 0.306^^^^ 


0.154lo-j;^^ < 0.225 

PI -,70+0:057 „ r,. 0+0. 089 
0.1'3_oo81 0.^43_o.204 
n OzL'S+0058 n onq+0.113 
u.zy:a_Q 097 "■'-'^''-0.280 


JW'500-7b.2-0.5r5oo 

-'W'l000-Tb.2_0.5r5oo 

A^2500-Tb.2-0.5r5oo 


0.327to;i36 < 0.503 0.273+^;^^^ 
0.267l:o:?l3 < 0.364 0.292lj;;^5° 
0.226to;?i7 < 0.295 0.328tJ!;^^^ 


0.288to:i33 < 0.280 
0.228+°-?^^ < 0.318 
0.199«;?»° < 0.411 


Msoo-Mgas 

Miooo-.'^gas 

-'^2500^A^gas 


0.123+^;^^^ < 0.260 < 0.226 
0.160+°°™ < 0.330 0.207+° ;;°^ 

n 9ZL1 +0.'077 „ .oc+0.190 n. „^o+o:i21 
U.^4i_Qn4 U.4db_o 410 U./Clci_o 230 


0.109+^-^15^ < 0.225 
0.147to:?o4 < 0-225 
0.2mto°l7 < 0.248 



Note. — Intrinsic scat ter, (TinMi in the mass-observable relations, using weak-lensing masses, for all 12 clusters and the 11 clusters 
excluding A1914 (Sec. 13.41 1. respectively. We refer to the mean of the posterior probability distribution of each parameter. Column (1): 
Scaling relations Columns (2-4): Intrinsic scatter for all 12 clusters, undisturbed and disturbed clusters, respectively. Columns (5-6): 
Intrinsic scatter for the 11 clusters and 6 disturbed clusters, excluding A1914, respectively. 



TABLE 4 
Multi-dimensional Fitting for covariance of intrinsic scatter of gas mass and temperature at A=500 



12 clusters 

(1) 



Flat Prior 

(A. 



U.Otiti 

0.086 

+0.055 

0.080 

> 0.185 
< 0.565 



Gaussian Prior 

^ 



'J-J^^J^-0.061 

206+°-°*3 
U.ZUD_o 051 

> 0.170 
< 0.597 



"3 
crt 



0.132 
0.213 



V{\r'\ > \r\) 



Note. — The model parameters in the covariance of intrinsic scatter for gas mass and temperature, derived from the multi-dimensional 
fitting for the full sample. We refer to the mean of the posterior distribution of each parameter. The lower limit is at a 68.3% confidence 
level. Column (1) : Model parameters. Column (2) : Results with the fiat prior Pprior = 1- Column (3) : Results with the Gaussian 
prior with the best-fit intrinsic scatter derived from the independent measurement for mass-observable relation. P(|r'| > \r\) denotes the 
probability that the correlation coefficient of the two random variables for twelve pair realizations is higher than the observed one. 



